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PREFACE 

This  Memorandum  reports  recent  results  of  continuing 
RAND  research  on  the  theoretical  and  computational  aspects 
of  ideal  chemical  systems  [1-13].  This  research  is  com¬ 
panion  to  studies  on  the  use  of  mathematical  models  to 
investigate  the  chemistry  of  physiological  subsystems 
[14-21].  Specifically,  the  Memorandum  deals  with  theo¬ 
retical  properties  of  the  type  of  ideal  chemical  system 
suggested  by  problems  involving  semipermeable  membranes 
in  biological  systems. 

The  discussions,  results,  and  informal  proofs  should 
be  intelligible  to  mathematically  literate  biochemists, 
physiologists,  research  physicians,  etc.  Some  of  the 
formal  proofs,  however,  require  more  mathematical  back¬ 
ground. 
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SUMMARY 


This  study  mathematically  investigates  equilibrium 
properties  of  certain  two-phase  chemical  systems.  The 
different  types  of  equilibrium  compositions  and  the  con¬ 
ditions  that  govern  them  are  determined.  The  effects  of 
linear  variation  in  the  chemical  inputs  are  studied;  e.g., 
the  effects  on  a  multisubstance  fluid,  separated  by  a 
semipermeable  membrane,  to  which  another  fluid  is  slowly 
added.  Among  other  results,  it  is  shown  that  the  rela¬ 
tive  sizes  of  the  two  phases  will  always  change  monoton- 
ically,  but  that  the  absolute  sizes  may  fluctuate. 
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SYMBOL  GLOSSARY 


bi(t)  =  bi  +  tAbi 


n 

b  =  I  b.  +  q  +  q 1 
i=l  1 

n 

b(t)  =  q (t)  +  q'(t)  +  F  b.(t) 

i=l 


ci 
c ! 

i 

Ab. 

1 


Aq' 


F 

g 

i 


Input  for  X^,  number  of 
moles  of  in  system. 

Total  number  of  moles  of 
Xt  in  S(0). 

Total  number  of  moles  of 

X.  in  S(t)  . 

1 

Total  number  of  moles  in 
system. 

Total  number  of  moles  in 
S(t). 

Gibbs  Parameter  for  X*. 

Gibbs  Parameter  for  X^1. 

Rate  at  which  number  of 
moles  of  X^  in  S(t)  is 
changed  per  unit  change  in  t. 

Rate  at  which  number  of  moles 
of  impermeable  material  in 
Phase  I  of  S(t)  is  changed 
per  unit  change  in  t. 

Rate  at  which  number  of  moles 
of  impermeable  material  in 
Phase  II  of  S(t)  is  changed 
per  unit  change  in  t. 

Free  energy. 

Characteristic  function. 

Index  of  permeable  species. 


Equilibrium  constant  for  the 
VI 

reaction  X.  ^  X,  . 

1  i 

Number  of  permeable  species. 

Number  of  moles  of  impermeable 
material  in  Phase  I. 

Number  of  moles  of  impermeable 
material  in  Phase  II. 

Number  of  moles  of  impermeable 
material  in  Phase  I  of  S(0). 

Number  of  moles  of  impermeable 
material  in  Phase  II  of  S(0). 

Number  of  moles  of  impermeable 
material  in  Phase  I  of  S(t). 

Number  of  moles  of  impermeable 
material  in  Phase  II  of  S(t). 

Partial  pressure  of  X„ 

Vapor  pressure  of  X^. 

Gas  constant. 

Variable  system  at  t. 

Variable  ranging  from  0  to  b. 

Root  of  the  characteristic 
function . 

Temperature . 

Used  in  connection  with  variable 
systems . 

Name  of  the  iun  permeable  species. 
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X1 

1 

X.  in  Phase  I. 

i 

x11 

1 

X.  in  Phase  II. 

l 

X. 

1 

Number  of  moles  of  X^  at 

l 

equilibrium. 

x! 

i 

Number  of  moles  of  X^  at 

l 

equilibrium. 

xt(t) 

Number  of  moles  of  X^  in 

l 

equilibrium  composition 

of  S (t) . 

x'(t) 

Number  of  moles  of  X^  in 

equilibrium  composition 

of  S(t) . 

xi 

x.  *  — 

l 

Mole  fraction  of  X*  at 

X 

equilibrium. 

x! 

* ,  i 

1  ’  X- 

Mole  fraction  of  X . ^  at 

l 

equilibrium. 

X.(t) 

X. (t)  - 

1  J(t) 

Mole  fraction  of  X*  in 

l 

equilibrium  composition 

of  S (t) . 

5!(t)  = 

1  i(t) 

T  J 

Mole  fraction  of  Xt"  in 

l 

equilibrium  composition 

of  S (t) . 

n 

x  =  I.  x.  +  q 
i=l  1 

Total  number  of  moles  in 

Phase  I  at  equilibrium. 

n 

x*  =  I  x!  +  q' 
i=l  1 

Total  number  of  moles  in 

Phase  II  at  equilibrium. 

xiv 


n 

x(t)  =  q ( t)  +  E  x.(t) 
1=1 


x'(t)  =  q'(t)  +  E  x! (t) 

1=1 


Total  number  of  moles  in 
Phase  I  in  equilibrium 
composition  of  S(t). 

Total  number  of  moles  in 
Phase  II  in  equilibrium 
composition  of  S(t). 
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I .  INTRODUCTION 


The  archetype  of  the  mathematical  system  explored  in 
this  Memorandum  is  a  multisubstance  fluid  partitioned  by 
a  single  semipermeable  membrane,  although  it  is  exemplified 
by  a  variety  of  other  chemical  equilibrium  problems  (e.g.  , 
a  gas  over  a  liquid).  More  explicitly,  the  Memorandum 
studies  two-phase,  constant-pressure-temperature,  chemical- 

* 

equilibrium  systems  that  satisfy  the  following  assumptions: 

1)  There  are  no  charged  species; 

2)  There  are  no  intra-phase  reactions--!. e. ,  the  only 
possible  "reactions"  are  the  migrations  of  molecules 
from  one  phase  to  another  phase; 

3)  The  substances  are  miscible  in  all  proportions. 

Although  such  simple  systems  rarely  occur  in  physiological 

applications,  their  study  is  important  for  several  reasons. 

First,  they  permit  a  more  thorough  mathematical  analysis 

than  more  complex  systems.  Second,  much  of  the  behavior 

of  the  latter  can,  at  least  qualitatively,  be  explained  in 

terms  of  the  behavior  of  the  simpler  systems--which,  from 

th 

one  point  of  view,  are  zero  -order  approximations  to  many 
physiological  systems.  Finally,  despite  their  apparent 

See  below  for  an  exact  definition  of  these  systems. 


-2- 


simplicity,  the  systems  we  study  display  a  wealth  of  com¬ 
plex  -rid  interesting  behavior. 

We  call  the  type  of  chemical  equilibrium  systems 
studied  in  this  Memorandum  primitive  chemical  equilibrium 
systems ,  or  more  briefly,  primitive  systems.  A  primitive 
system  will  involve  two  types  of  chemical  substances:  im¬ 
permeable  substances,  whose  molecules  cannot  pass  through 
the  membrane,  and  permeable  substances,  whose  molecules 
can  pass  through  the  membrane. 

Note  that  for  primitive  systems  the  only  practical 
effect  of  the  presence  of  impermeable  substances  is  on  the 
mole  fractions  of  the  permeable  substances.  For  this 
reason,  there  is  no  loss  of  generality  in  aggregating  all 
the  impermeable  substances  to  a  single  impermeable  sub¬ 
stance.  (This  would  not  necessarily  be  the  case  if  assump¬ 
tions  1  and  2  above  did  not  hold.)  We  shall  deal  therefore 
with  a  system  containing  n  permeable  substances  ,  named 
, . . . ,X^ , . . .  ,X  ,  and  one  unnamed  impermeable  substance. 

And  since  we  can  always  say  that  the  amount  of  impermeable 
material  present  happens  to  be  zero,  the  fact  that  we  have 
specified  an  impermeable  substance  does  not  preclude  us  from 
studying  problems  in  which  there  are  no  impermeable 


substances . 
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We  designate  the  two  fluids  separated  by  the  membrane 
as  Phase  I  and  Phase  II.  Note  that  one  phase  might  lie 
above  or  to  the  side  of  the  other  phase;  or  the  membrane 
might  be  a  surface  enclosing  a  volume,  the  two  phases 
being  the  inside  and  outside  of  the  membrane;  or  the 
membrane  might  be  several  such  closed  surfaces  and  the 
two  phases  might  be  the  collective  interiors  and  the 
exterior.  An  example  is  the  blood  plasma  and  the  interiors 
of  the  red  blood  cells. 

When  the  substance  X.  occurs  in  Phase  I,  we  shall  call 

L 

I  II 

it  X.;  when  it  occurs  in  Phase  II  we  shall  call  it  X.  . 

1  1 

Thus,  the  migration  of  X^  across  the  membrane  is  repre¬ 
sented  by  the  "chemical  reaction" 

X1  *  X11  .  (1) 

li 

Let  x.  and  x!  be  the  number  of  moles  of  X'f  and  X* 

l  i  -  i  l 

respectively;  let  q  and  q'  be  the  number  of  moles  of  im¬ 
permeable  material  in  Phase  I  and  Phase  II,  respectively; 
and  let  x  and  x'  be  the  total  number  of  moles  in  Phase  I 
and  Phase  II,  respectively.  That  is, 

n 

=  q  + 1  \  ; 
i=l 


x 


(2) 


up¬ 


l¬ 


and 


n 


x'  =  q'  +  £  x[  . 

i=l 


(2') 


Let  x,  and  x!  be  the  mole  fractions  of  and  X^, 
i  l  -  i  i 

respectively.  That  is, 


x^  =  —  for  i=l , .  .  .  ,  n  ; 
x 


(3) 


and 


.  i  _ 


x! 

l 
'  i 


for  i=l, . . . ,n  . 


(3') 


Let  be  the  total  number  of  moles  of  X^  in  the  system. 


That  is, 


x.  +  x!  =  b. 
ill 


for  i=l , .  .  .  ,n 


(4) 


We  call  b^  the  input  for  X^.  Let  b  be  the  total  number  of 
input  moles.  That  is, 


If  x  =  0  then  x.  is  undefined,  and  if  x'  =0  then 

l  ’ 

xj^  is  undefined. 
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n 

b  =  q  +  q'  +  £  .  (5) 

i=l 

Equations  (2),  (2'),  (4),  and  (5)  show  that 

x  -f  x '  =  b  .  (6) 

We  assume  that  the  system  is  in  equilibrium  when  its 

free  energy,  as  defined  below,  is  minimized  with  respect 

to  the  variables  x,,...,x  ,x',...,x’,  constrained  by  Eqs. 

1’  n’l  n  j 

(4)  above  and  the  following: 

x^  ^  0  for  i=l,...,n  ;  (7) 


and 


x!  ^  0  for  i=l , . .  .  ,  n 
l  5  ’ 


(7') 


Finally,  we  assume  that  the  free  energy  is  given  by 


an  equation  of  the  form: 

n 

F  = 

i=l 


=  ][  (c”i  +  log  x^)+  q  log  ^ 


(8) 


+  I  Xi  (^i  +  108  XD  +  q'  log 

i=l  V  ^  x' 


-6- 


Here  c^  is  the  Gibbs  parameter  for  X^,  and  is  the  Gibbs 
parameter  for  x|P,  and  x^,  x|,  x,  and  x*  are  defined  in 
terms  of  the  and  x|  by  means  of  Eqs.  (2),  (21),  (3), 
and  (3'). 

Several  comments  concerning  the  equation  for  the  free 
energy  function  are  in  order: 

1)  The  right-hand  side  is  equal  to  what  is  classically 
called  the  Gibbs  Free  Energy  divided  by  the  constant 
factor  RT,  where  R  is  the  gas  constant  and  T  is 

the  absolute  temperature. 

2)  Its  form  assumes  that  the  system  is  "ideal"--to 

the  extent  that  the  so-called  mole  fraction  activity 
coefficients  are  independent  of  composition.  While 
not  true  in  practice,  this  assumption  is  often  suf¬ 
ficiently  accurate  over  a  wide  range  of  compositions. 

3)  At  first  glance,  the  parameters  c^  and  c!,  would  seem 
to  be  necessarily  the  same.  However,  on  further 
reflection,  this  is  seen  to  be  true  only  when  the 
membrane  is  chemically  symmetric  and  inert,  when 
the  pressures  and  temperatures  on  both  sides  are 
equal,  and  when  the  compositions  of  the  two  phases 
are  sufficiently  similar  that  the  mole  fraction 
activity  coefficients  are  equal.  None  of  these  as¬ 
sumptions  need  hold.  Furthermore,  if  we  take  as 
our  model  of  a  primitive  system  a  liquid-gas  system 
(instead  of  two  liquids  separated  by  a  semipermeable 
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membrane) ,  then  there  is  absolutely  no  reason  to 
suppose  that  =  c|. 

4)  Equation  (8)  as  written  defines  the  free  energy 
only  if  each  x^  >  0,  each  x^  >  0,  and  q  >  0  and 
q1  >  0.  The  free  energy  must  also  be  defined 
when  one  or  more  of  these  quantities  are  zero. 

This  can  be  done;  for  a  discussion  of  this  and 
related  matters,  see  Ref.  8,  Theorem  8.3. 

It  turns  out  [11]  that  the  behavior  of  a  primitive 

system  depends  on  the  Gibbs  parameters  only  to  the  extent 

that  it  depends  on  the  differences,  c.  -  cl.  Thus,  if 

c .  -  c  ! 

we  define  =  e  ,  then  a  primitive  system  can  be 

specified  by  the  number  of  permeable  substances,  n;  the 
inputs,  , . . . ,bn,q ,q' ;  and  the  constants ,  k^,...,kn. 

For  equilibrium  compositions  in  which  x  >  0  and  x'  >  0, 

•k 

we  must  have: 


• 

—  =  k .  for  1*1, .  .  .  ,n  .  (9) 

A  L 

X. 

1 


Equation  (9)  expresses  the  law  of  mass  action  for  the 


"reaction"  (1).  This  motivates  us  to  refer  to  the  k^  as 


See  Ref.  8,  p.  44,  Theorem  10.3. 
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equiltbrium  constants.  (The  k^  are  sometimes  called 
’’partition  coefficients".)  The  equilibrium  constants  are 
more  convenient  parameters  for  discussing  the  behavior  of 
chemical  equilibrium  systems  than  are  the  Gibbs  parameters. 

Conversely,  any  composition  xi  > • • • >xn >9 >  x^,...,x^,q' 
which  satisfies  x  >  0,  x'  >0,  the  law  of  mass  balance 
Eq.  (4),  and  the  law  of  mass  action  Eq.  (9)  is  an  equilibrium 

•k 

composition  of  the  system  in  question. 

Note  that  0  <  k^  <  «.  From  Eq.  (9),  it  would  appear 
that  if  we  allowed  k^  =  0,  then  we  might  regard  X^  as  being 
impermeable  material  in  Phase  I;  and  that  if  we  allowed 
k^  =  00 ,  then  we  might  regard  as  being  impermeable  ma¬ 
terial  in  Phase  II.  The  precise  sense  in  which  this  is 
true  is  treated  for  general  chemical  equilibrium  systems 
in  Ref.  11.  In  this  Memorandum,  we  continue  to  assume 
0  <  k^  <  00 ,  and  to  treat  impermeable  species  in  a  special 
way. 

•k 

Note  that  since  we  have  defined  equilibrium  in  terms 
of  the  free  energy  function,  this  statement  requires  proof. 
Also,  note  that  it  says  nothing  about  equilibrium  composi¬ 
tions  for  which  x  =  0  or  x'  =0.  See  Ref.  8,  p.  44,  Theorem 
10.3,  for  a  proof. 
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II,  THE  CHARACTERISTIC  FUNCTION 

Assume  that  x^,...,xn>q,  x|,...,x^,q'  form  an  equi¬ 
librium  composition  and  that  x  >  0  and  x'  >  0.  If  in  Eq . 
(9)  we  replace  x^  and  x^  by  their  values  from  Eqs.  (3) 
and  (3')  and  apply  Eq.  (4),  then 

b .  x 

x  =  :--1  J  (10) 

x  +  k .  x' 

l 


and 


k  ,b .  x' 

L  1 

x  +  k.x' 

L 


(10') 


If  we  replace  x^  in  Eq.  (2)  by  its  value  from  Eq.  (10), 
then 


.  +  7  t — Hr  - 1  •  (u) 

x  .  ,  x  +  k.x' 
i=l  l 


If  we  use  Eq.  (6),  then  Eq.  (11)  can  be  rewritten  as  g(x)  =  0 
where 


n 


g(s) 


i=l 


b. 

l 


s  +  k^(b  -  s) 


1  . 


(12) 
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We  call  g  the  characteristic  function  of  the  system. 

We  have  shown  that,  if  x  >  0  and  x'  >0,  and  if  the 

and  x^  define  an  equilibrium  composition,  then  x  is  a 
root  of  the  characteristic  function.  The  characteristic 
function  plays  an  important  role  in  the  theory  of  primi¬ 
tive  systems. 

Now  drop  the  assumption  that  x  >  0  and  x*  >  0.  As¬ 
sume  first  that  x  =  0.  Then  each  x^  =  0,  and  hence,  by 
Eq.  (4),  each  x|  =  b^.  Inspection  of  Eqs.  (10)  and  (10') 
shows  that  they  are  satisfied;  similarly,  they  are  satis¬ 
fied  if  x1  =  0.  It  follows  as  above  from  Eqs.  (10)  and 
(10')  that  g(x)  =  0  for  these  cases.  Thus,  we  have  shown 
that  for  any  primitive  chemical  equilibrium  system  Eqs, 

(10)  and  (10*)  are  satisfied  and  x  is  a  root  of  the  charac¬ 


teristic  function. 
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III.  TYPES  OF  EQUILIBRIUM  COMPOSITION 

This  section  states  and  proves  our  principal  results 
concerning  the  types  of  equilibrium  compositions  possible 
in  a  primitive  chemical  equilibrium  system,  the  conditions 
for  each  type  to  occur,  the  properties  of  the  character¬ 
istic  function  of  a  primitive  system,  and  the  relationships 
between  the  characteristic  function  and  the  equilibrium 
compositions.  This  requires  the  definition  of  two 
conditions . 

"  b. 

Condition  I:  Either  q>0  or  V  >  b  . 

1-1 1 

n 

Condition  II:  Either  q'  >  0  or  £  kJt^  >  ^  • 

i=l 

Observe  that  if  Condition  I  is  false,  then  we  must  have 
n  b . 

q  =  0  and  I  s;  b.  Similarly,  if  Condition  II  is  not 
i=l  i  n 

satisfied,  then  q*  =  0  and  E  k.b.  £  b. 

i=l 

THEOREM  1 

Consider  a  primitive  chemical  equilibrium  system 
with  n  permeable  species  having  inputs  b^,...,b  and 
equilibrium  constants  k^,...,kn>  and  with  an  impermeable 
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species  having  input  q  in  Phase  I  and  input  q'  in 

n 

Phase  II.  Let  b  =  E  b.  +  q  +  q1. 

i=l  1 

(i)  The  characteristic  function  g(s)  of  the  system 

*  _ 

is  defined  and  convex  on  the  interval  0  <  s  <  b. 

(ii)  If  both  Condition  I  and  Condition  II  are 
satisfied,  then  g  has  exactly  one  root,  Sq,  in  that 
interval,  and  the  system  has  exactly  one  equilibrium 
composition,  given  by 


x. 

1 


biso 


so  +  ki(b 


for  i=l , . .  .  ,  n  ; 


and 


(14) 


kibi(b  -  V 


+  k. (b  - 

l 


for  i=l 


,n 


(15) 


Furthermore , 


x 


(16) 


A  function  is  convex  if  its  graph  is  such  that  it 
lies  below  all  its  chords.  For  example,  a  function  having 
an  everywhere-positive  second  derivative  is  convex.  More 
detailed  discussions  of  convexity  appear  in  Ref.  23,  pp. 
149-156. 
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x' 


(17) 


A 

X. 

1 


b. 

L 


+  k. (b  - 

l ' 


for 


i=l 


,n 


(18) 


and 


A 


k.b. 

1  L 


s0  +  k.(b  -  s0) 


for 


i=l , . . . , n  . 


(19) 


(iil)  If  either  Condition  I  or  Condition  II  is 
satisfied,  but  not  both,  then  g  has  no  roots  in  the 
interval  0  <  s  <  b. 

(iv)  If  Condition  I  is  satisfied  but  not  Condition  II, 
then  the  system  has  exactly  one  equilibrium  composition, 
given  by  x^  =  b^  and  xl^  =  0,  for  i=l,...,n. 

(v)  If  Condition  II  is  satisfied  but  not  Condition  I, 
then  the  system  has  exactly  one  equilibrium  composition, 
given  by  x^  =  0  and  x^  =  b^,  for  i=l,...,n. 

(vi)  If  neither  Condition  I  nor  Condition  II  is 
satisfied,  then  q  =  q'  =0,  and  b^  =  0  for  all  i, 
i=l,...,n,  for  which  k^  ^  1.  Conversely,  if  q  =  q'  =0 
and  b^  =  0  for  all  i,  i=l,...,n,  for  which  k^  ^  1,  then 
neither  Condition  I  nor  Condition  II  is  satisfied. 
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In  this  case,  the  system,  with  trivial  exceptions, 
has  an  infinite  number  of  equilibrium  compositions, 
given  by  =  tb^  and  xl^  =  (1  -  t)b^,  for  i=l,...,n, 
where  each  t,  O^t^l,  gives  a  different  equilibrium 
composition. 

Proof  of  Theorem  1; 

If  any  b^  =  0,  then  x^  =  x|  =  0  for  every  equilibrium 
composition.  Thus,  if  we  ignore  every  species  for  which 
b^  =  0  and  appropriately  decrease  n  (the  number  of  perme¬ 
able  species) ,  then  the  characteristic  function  is  un¬ 
changed,  the  satisfaction  of  Conditions  I  and  II  is  un¬ 
changed,  and  none  of  the  statements  of  the  theorem  are 
affected.  We  may  accordingly  assume  that  each  b^  >  0. 
Also,  the  theorem  is  trivially  true  if  n  =  0,  so  we  will 
assume  that  n  >  0. 

The  characteristic  function  as  defined  by  Eq.  (12) 
is  obviously  defined  on  0  <  s  ^  b  if  q  >  0,  and  on 
0  ^  s  s  b  if  q  =  0.  Its  first  two  derivatives  are  easily 
calculated : 


If  every  b^  =  0,  or  if  n  =  0,  then  there  is  a  unique 
equilibrium  composition. 


-15- 
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=  .  £_  .  V 

ds  2  L 


b . (1  -  k.) 

L  1 


L  r 

i=l  s(l  -  k.)  +  k.b 

L  i  i  J 


(20) 


-  la 
,  2  3 

ds  s 


n 

♦’  l 


b  .  (1  -  k.) 

i  l 


L  r 

1=1  s ( 1  -  k.)  +  k.b 

L  1  1  J 


(21) 


The  consequences  will  be  detailed  in  several  Lemmas.  First, 
it  is  apparent  from  Eq .  (21)  and  from  the  definition  of 
g,  Eq.  (12),  that 
Lemma  1 

The  characteristic  function,  g,  is  continuous 
and  convex  on  0  <  s  b.  It  is  strictly  convex  every¬ 
where  on  this  interval  unless  each  k.  =  1  and  q  =  0, 

l 

in  which  case  it  is  constant  and  equal  to  -  . 

b 

By  examining  the  behavior  of  g  near  0,  it  follows  that 
Lemma  2 

If  Condition  I  is  satisfied,  then  g(c)  >  0  for 
all  sufficiently  small  positive  c.  If  Condition  I  is 
not  satisfied,  then  g  is  defined  at  0,  and  g(0)  £  0. 

By  examining  the  behavior  of  g  and  its  derivative  at  b 


it  follows  that 
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Lemma  3 

If  Condition  II  is  satisfied,  then  g(b  -  e)  <  0 

for  all  sufficiently  small  positive  e.  If  Condition  II 

is  not  satisfied,  then  g(b)  =  0  and  ^  £  0. 

s=b 

It  then  follows  from  Lemmas  1,  2,  and  3,  that 
Lemma  4 

If  Conditions  I  and  II  both  hold,  then  g  has 
exactly  one  root  iri  the  interval  0  s  <  b;  if  exactly 
one  of  these  two  conditions  hold,  then  g  has  no  roots 
in  that  interval;  if  neither  condition  holds,  then  g 
vanishes  identically  and,  hence,  every  =  1. 

Lemma  5 

The  equilibrium  solutions  of  the  system  in  which 
x  >  0  and  x'  >  0  are  in  one-to-one  correspondence  with 
the  roots  of  g,  in  0  <  s  <  b.  This  correspondence  is 
defined  by  Eqs .  (14)  and  (15)  and  it  satisfies  Eqs . 

(16),  (17),  (18),  and  (19). 

To  prove  Lemma  5,  the  discussion  of  Sec.  II  shows  that 
every  equilibrium  composition  with  x,  x'  >  0  satisfies 
Eqs.  (14)  through  (19)  for  some  root,  Sq.  Conversely, 
given  any  root,  Sq,  we  may  define  x^  and  x|  by  means  of 
Eqs.  (14)  and  (15).  Equations  (16),  (17),  (18),  and  (19) 
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will  then  be  satisfied,  and  also  x.  >  0,  x!  >  0,  x.  + 

ill 

A  A 

x!  =  b.,  and  x!  =  k.x..  Thus,  the  x.  and  x!  are  positive 
li  ill  l  l  r 

and  satisfy  the  laws  of  mass  balance  and  of  mass  action 
for  the  system.  Hence,  the  x.  and  x!  as  defined  from  sn 

J  ’ii  0 

constitute  an  equilibrium  composition.  Since  Sq  =  x 
(Eq.  (16)),  the  correspondence  is  one-to-one. 

Lemma  6 

If  Conditions  I  and  II  are  both  satisfied,  then 
the  only  equilibrium  composition  is  the  one  defined 
by  Lemma  5,  associated  with  the  unique  (Lemma  4)  root 
of  g  in  0  <  s  <  b. 

To  prove  Lemma  6,  we  need  only  show  that  the  system 
can  have  no  equilibrium  composition  with  x  =  0  or  x'  =0. 

But  the  only  equilibrium  composition  with  x  =  0  is  given 
by  x^  =  0  and  x|  =  b^,  and  the  only  equilibrium  composition 
in  which  x*  =  0  is  that  for  which  x!  =  0  and  x.  =  b. .  But 

l  ii 

if  the  system  had  more  than  one  equilibrium  composition, 
it  would  have  infinitely  many  equilibrium  compositions  and, 
hence,  infinitely  many  equilibrium  compositions  with  x  >  0 
and  x'  >0,  contradicting  Lemmas  4  and  5. 


* 

See  Ref.  8,  p.  38,  Lemma  9.3. 
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Lemma  7 

If  Condition  I  is  satisfied,  but  not  Condition  II, 
then  there  is  exactly  one  equilibrium  composition.  It 
is  given  by  x.  =  b.  and  x!  =  0. 

b  J  l  l  l 

To  prove  this,  observe  that  by  Lemmas  4  and  5  every 
equilibrium  composition  must  have  x  =  0  or  x'  =0.  That 
is,  the  only  possible  equilibrium  compositions  are  those 
given  by  x.  =  0  and  x|  =  b^,  or  x^  =  b^  and  x|  =  0.  Exactly 
one  of  these  two  compositions  must  be  an  equilibrium  com¬ 
position,  because  the  system  possesses  at  least  one  equi- 

-k 

librium  composition,  and  if  there  were  two  there  would 

have  to  be  infinitely  many,  as  in  the  proof  of  Lemma  6. 

n 

If  q  >  0  then,  since  x  =  q  +  I  b.  >  0,  the  equilibrium 

i=l  1 

composition  must  be  given  by  x^  =  b^  and  x|  =  0,  as  asserted 
in  the  lemma.  To  show  that  this  is  also  the  case  when 
q  =  0,  let  us  modify  the  system,  producing  a  new  system 
J (q)  by  changing  q.  Thus  the  original  system  is  J(0). 

Since  ^f(0)  satisfies  Condition  I  but  not  II,  so  must  A(s\) 
satisfy  Condition  I  but  not  II  for  q  >  0.  But,  for  q  >  0, 
we  have  already  shown  that  ^f(q)  has  the  unique  equilibrium 
composition  x^(q)  =  b^,  x|(q)  =  0.  Hence 

See  Ref.  8,  p.  42,  Theorem  9.9. 
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11m  x  (q)  =  b  i=l , . . . , n 

q-0 

and 

lim  x! (q)  =  0  1=1,..., n  . 

q-0+  1 

By  the  continuity  results  of  Ref.  11,  it  follows  that 
Xi  =  ^i  anc^  x[  =  0  define  the  equilibrium  composition  for 
^  (0) ,  as  required. 

Lemma  8 

If  Condition  II  but  not  Condition  I  is  satisfied, 
then  the  system  has  exactly  one  equilibrium  composition, 
given  by  x|  =  0  and  =  b^. 

The  proof  of  Lemma  8  is  similar  to  the  proof  of  Lemma  7, 
or  follows  from  Lemma  7  by  reversing  the  roles  of  Phases  I 
and  II  and  replacing  the  k^  by  their  reciprocals. 

All  that  remains  to  complete  the  proof  of  Theorem  1 
is  to  observe  that  if  q  =  q'  =0  and  every  =  1  then 
clearly  neither  Condition  I  nor  Condition  II  can  hold,  and 
to  discuss  the  nature  of  the  solutions  in  this  case.  By 
Lemmas  4  and  5,  the  compositions  described  in  the  theorem 
for  this  case  are  equilibrium  compositions  for  0  <  t  <  1. 


-20- 


But  the  equilibrium  compositions  are  known  to  form  a  closed 
* 

set.  Thus,  the  compositions  described  are  equilibrium 
compositions  for  0  *  t  £  1 .  That  there  are  no  other 
equilibrium  compositions  also  follows  from  Theorem  9.2 
in  Ref.  8. 

See  Ref.  8,  p.  38,  Theorem  9.2. 
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IV.  NUMERICAL  COMPUTATION 

The  practical,  numerical  computation  of  the  equilibrium 
compositions  of  a  primitive  system,  given  its  equilibrium 
constants  and  inputs,  can  be  accomplished  either  by  means 
of  digital  computer  programs  designed  for  the  solution  of 
general  chemical  equilibrium  programs  (e.g.,  Refs.  1,  5, 

10),  or  with  the  special  techniques  described  in  this 
section . 

A  procedure  for  the  determination  of  the  equilibrium 
composition  of  a  primitive  system  would  first  determine  the 
values  of  the  and  the  x|.  From  these,  other  quantities 

A  A 

such  as  the  mole  fractions,  x. ,  x! ;  volumes;  concentrations 
in  other  scales,  etc.,  could  be  calculated.  The  x.  and  the 

l 

x|  may  be  calculated  from  x  and  x'  through  Eqs.  (10)  and 
(10')  of  Sec.  II.  Since  x'  =  b  -  x,  the  calculation  re¬ 
duces  to  a  technique  for  calculating  x.  To  do  this,  one  can 
evaluate  the  validity  of  Condition  I  and  the  validity  of 
Condition  II  of  Sec.  Ill,  and  then  use  Theorem  1  to  deter¬ 
mine  that  x  =  0,  or  x  =  b,  or  x  is  an  indeterminate  be¬ 
tween  0  and  b,  or  that  x  is  the  unique  root  of  the  charac¬ 
teristic  function  in  the  interval  0  <  s  <  b.  The  charac¬ 
teristic  function  being  convex  on  0  <  s  <  b  guarantees  that 
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Newton's  method  will  converge  to  the  root  as  long  as  we 
begin  with  an  initial  value  in  the  interval  0  <  s  <  b. 

Although  g(s)  =  0  can  be  rewritten  as  a  polynomial 
equation,  this  is  of  little  practical  use  except  for  very 
small  systems.  For  example,  the  primitive  system  treated 
in  Ref.  21  has  a  characteristic  function  whose  root  can  be 
determined  by  solving  a  quadratic  equation,  and  thus  can 
be  solved  analytically. 

i< 

See,  for  example,  Ref.  22,  pp.  81-82,  or  any  good 
numerical  analysis  textbook. 
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V,  GAS-LIQUID  SYSTEMS 


In  this  section,  to  illustrate  the  application  of  our 
results  to  other  than  membrane  phenomena,  we  consider  a 
gas-liquid  system  of  n  substances,  X^,...,X  ,  at  a  pressure 
P  and  temperature  T. 

We  assume  that  under  these  conditions  each  substance 
will  be  distributed  between  a  single  liquid  phase  and  a 
single  gas  phase.  Thus,  we  are  assuming  that  the  substances, 
in  both  liquid  and  solid  form,  are  completely  miscible  in  all 
proportions,  and  that  no  solid  precipitates  are  formed. 

We  also  assume  that  the  substances  are  involved  in  no 
chemical  reactions  other  than  evaporation  and  condensation. 

We  shall  show  that  if  this  type  of  system  satisfies  certain 
further  assumptions  it  is  a  primitive  chemical  equilibrium 
system.  Thus  all  our  results  concerning  primitive  systems 
will  be  applicable  to  this  type  of  gas-liquid  system. 

We  designate  the  liquid  phase  Phase  I,  and  the  gas 
phase  Phase  II,  and  let  x^,  xl^,  iL,  x^ ,  b^,  x,  x' ,  etc., 
have  their  usual  meanings. 

Let  the  vapor  pressure  of  at  the  temperature  T  be 
p?.  We  assume  that  each  p°  is  positive  and  finite;  i.e., 
we  assume  that  each  X^  is  neither  completely  volatile  nor 
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completely  non-volatile.  This  is  analogous  to  the  con¬ 
dition,  for  a  membrane  system,  that  there  are  no  impermeable 
substances.  We  assume  that  the  liquid  obeys  Raoult's  law-- 
which  we  take  to  be  that  the  partial  pressure,  p^,  of 
is  equal  to  the  product  of  the  mole  fraction,  x^,  of  X* 
and  the  vapor  pressure,  p? ,  of  X^: 


(30) 


We  assume  that  the  gas  obeys  Dalton's  law  of  partial  pres- 
sures--which  we  take  to  be  that  the  partial  pressure,  p^, 
of  X*1  is  equal  to  the  product  of  the  mole  fraction,  x|, 
of  xjf^  and  the  total  pressure,  P: 

=  Px[  .  (31) 


Since  at  equilibrium  X^  must  have  the  same  partial  pressure 
in  liquid  as  in  gas,  we  have 


P 


x . 

L 


(32) 
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Thus  ,  the  system  is  a  primitive  system  with  the  equilibrium 
constant,  k^,  of  given  by 


o 


As  an  example  of  the  application  of  our  general  re¬ 
sults  to  gas-liquid  systems,  let  us  find  the  conditions 
for  a  non-unique  solution.  According  to  Theorem  1  of 
Sec.  Ill,  there  will  be  a  non-unique  equilibrium  composi¬ 
tion  if  and  only  if  =  1  for  every  i  for  which  b^  £  0, 

or,  applying  Eq.  (33),  if  and  only  if  p?  =  P  for  every  i 
for  which  b^  £  0.  In  other  words,  there  will  be  a  non¬ 
unique  equilibrium  composition  at  pressure  P  and  tempera¬ 
ture  T  precisely  if,  for  every  substance  present  in  any 
positive  amount,  the  vapor  pressure  is  P  when  the  tempera¬ 
ture  is  T;  i.e. ,  precisely  if  the  boiling  point  of  every 
substance  present,  at  the  pressure  P,  is  T. 

As  another  example,  let  us  find  the  circumst  inces 
under  which  a  liquid  phase  can  exist  at  equilibrium  without 
the  formation  of  a  gas  phase;  i.e.,  the  circumstances  re¬ 
quired  for  either  a  non-unique  solution  or  a  unique  solution 
with  no  gas  phase  present.  Applying  Theorem  1  of  Sec.  Ill 
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shows  that  these  are  precisely  the  circumstances  under 

bi 

which  Condition  II  fails--that  is  (because  x.  =  — ): 

1  e 


u 


l  kixi  ‘  1  : 

i-1 


or,  replacing  by  its  value  from  Eq .  (33)  and  using 
Eq.  (30), 

n 

£  Pi  -  P  •  (34) 

i-1 


Thus,  we  have  obtained  the  familiar  condition  that 
a  liquid  phase  can  exist  at  equilibrium  without  the  forma¬ 
tion  of  a  gas  phase  if  and  only  if  the  sum  of  the  partial 
pressures  of  its  constituents  does  not  exceed  the  hydro¬ 
static  pressure. 

Similarly,  we  can  obtain  the  condition  for  a  gas 
phase  to  exist  at  equilibrium  without  the  formation  of  a 
liquid  phase: 


n 
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VI.  INTRODUCTION  TO  VARIABLE  SYSTEMS 

Previous  sections  have  investigated  the  properties  of 
a  fixed  primitive  chemical  equilibrium  system.  This  and 
the  following  sections  study  the  manner  in  which  the  equi¬ 
librium  properties  of  a  variable  system  vary  when  the  in¬ 
puts  are  varied.  By  a  variable  system  we  mean  specifically 
a  primitive  system  with  n  permeable  species  and  fixed  equi¬ 
librium  constants,  but  with  inputs  that  vary  linearly  with 
a  parameter  t  which  may  or  may  not  denote  time.  In  any 
case,  we  emphasize  that  we  are  not  studying  kinetic  effects-- 
at  each  value  of  t,  the  system  is  allowed  to  come  to 
equilibrium. 

Thus,  the  amount  of  input  is  given  by  a  function 

b . (t)  =  b?  +  t£b .  . 

L  l  1 

Similarly,  the  amount  of  impermeable  material  in  Phase  I  is 

q ( t )  =  q°  +  t£q  , 

and  the  amount  of  impermeable  material  in  Phase  II  is 


q'(t)  -  q°’  +  tA,1  . 
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We  assume  that  Ab^  0,  Aq  *  0,  Aq  '  --  0,  and  t  >  0. 
Although  some  of  our  results  are  valid  without  this  as¬ 
sumption,  it  is  more  convenient  to  make  it  throughout. 

Thus,  each  t  -  0  presents  a  primitive  chemical  equi¬ 
librium  system,  which  we  shall  call  S(t).  If  Aq  =  Aq '  =  0 
n 

and  I  Ab.  =  1,  then  S(t)  may  be  interpreted  as  the  result 
i=l  1 

of  adding  t  moles  of  a  mixture  of  the  permeable  substances 

n 

with  mole  fractions  Ab.  , .  . . ,Ab  to  S(0).  If  I  Ab.  £  1, 

1  n  .  .  l 

i=l 

then  the  same  type  of  interpretation  holds,  except  that  the 
Ab^  are  no  longer  mole  fractions  and  t  is  no  longer  measured 
in  moles.  Finally,  if  Aq  £  0  or  Aq'  ^  0 ,  we  may  think  of 
adding  a  mixtrre  of  both  impermeable  and  permeable  sub¬ 
stances  to  S(0),  taking  care  to  distinguish  between  im¬ 
permeable  material  destined  for  Phase  I  and  impermeable 
material  destined  for  Phase  II. 

This  section  concludes  with  a  sequence  of  examples  of 
variable  primitive  chemical  equilibrium  systems  (leaving 
general  results  to  later  sections) . 

Example  1 

Consider  an  inert  symmetric  membrane  permeable  to 
l^O  but  not  to  some  macromolecule,  and  capable  of  main¬ 
taining  a  hydrostatic  pressure  gradient.  Suppose  we  place 
one  mole  of  the  macromolecule  on  one  side  (Phase  I)  of  the 
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membrane,  fix  the  hydrostatic  pressure  of  both  phases  so 
that  the  pressure  of  Phase  I  is  higher  than  the  pressure 
of  Phase  II,  i  j1  oth  phases  at  the  same  temperature,  and 
put  t  moles  of  H^O  into  the  system.  Thus,  we  have  q°  =  1, 
Aq  =  q°'  =  Aq '  =  0,  n  =  1,  Xx  =  H20,  b°  =  0,  &>l  =  1,  and 
>  1  (because  of  the  hydrostatic  pressure  gradient) .  We 
see  that  S(t)  satisfies  Condition  I  of  Sec.  Ill  for  all 
t  2  o,  and  satisfies  Condition  II  if  and  only  if  t  >  -j^y 
moles . 


It  follows  from  Theorem  1  that  S(t)  has  a  unique  equi¬ 
librium  composition  for  all  t  ^  0,  and  that  for  0  ^  t  £  j—y, 
x^  =  0  (i.e.  ,  all  the  1^0  will  be  drawn  to  the  side  of  the 
membrane  having  the  macromolecule) .  But  for  t  >  r~~y , 

V1 

x^  =  ~—y  and  x|  =  t  -  ■k-1-y;  i.e.  ,  after  y^-y  m°les  of 

H^O  have  been  added,  all  additional  H^O  will  be  drawn  to 
the  side  not  containing  the  macromolecule. 
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but  S(t)  has  t  moles  of  X^.  We  will  first  describe 
several  realizations  of  such  a  system: 

1)  Consider  a  gas-liquid  system  of  the  type  de¬ 
scribed  in  Sec.  V,  with  two  substances  and 
X^ ,  each  at  the  temperature  T.  The  vapor  pres¬ 
sure  of  X^  exceeds  that  of  X^,  and  the  total 
pressure,  P,  of  the  system  is  held  at  a  fixed 
value  between  the  vapor  pressure  of  X^  and  the 
vapor  pressure  of  X^. 

2)  Consider  two  substances,  X^  and  X^ ,  both  perme¬ 
able  to  an  inert  symmetric  membrane  which  is 
both  rigid  (i.e.,  capable  of  supporting  a  pres¬ 
sure  difference)  and  a  heat  insulator  (i.e., 

/V 

capable  of  supporting  a  temperature  difference). 
Since  and  depend  in  different  ways  on  the 
temperatures  and  pressures,  it  is  quite  possible 
that  k^  >  1  >  k^. 

3)  Consider  a  metabolical ly  active,  non-symmetr ic 
membrare  which  "pumps”  X^  (which  might  be  l^O) , 

/V 

Such  membranes  can  be  obtained  from  the  same  scientific 
supply  houses  which  furnish  frictionless  pistons  and  com¬ 
pletely  ionizable  electrolytes. 
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from  Phase  I  to  Phase  II  and  which  "pumps" 

(which  might  be  glucose  )  from  Phase  II  to 
Phase  I. 

We  can  apply  Theorem  1  to  this  variable  system  to  see 
that,  for  all  t  ^  0,  it  will  have  a  unique  solution.  For 


Phase  II  will  be  vacuous.  For 


1-k  k  l-k9 

■■■  <  t  <  — —  —  ■ '  ■ 

V1  k2  V1  * 

neither  phase  will  be  vacuous.  For 


t  £ 


kL  l-k2 

k2  kl_1 


Phase  I  will  be  vacuous.  Figure  1  shows  a  plot  of  the 
composition  of  the  system  versus  t,  for  the  case  k^  = 

k 

2  4' 

•k 

Note  that  we  must  then  consider  only  those  values 
of  t  for  which  no  glucose  precipitate  will  be  formed. 


CO|CM 


(phase  II  empty)  (neither  phase  empty)  (phase  T  empty) 

Fig.l — Variation  of  composition  in  the  system  of  Example  2 
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Example  3 

Consider  the  result  of  equilibrating  .04  moles  of 
methane  and  .96  moles  of  normal  octane  with  t  moles  of 
a  mixture  of  normal  hexane  (with  a  mole  fraction  of  .96) 
and  normal  eicosane  (with  a  mole  fraction  of  .04)  at  a 
temperature  of  40°C  and  a  pressure  of  120  mm  of  Hg,  as¬ 
suming  that  Raoult's  law  and  Dalton's  law  of  partial 
pressures  are  satisfied.  We  have  n  =  4,  =  methane, 

X 2  =  normal  hexane,  X^  =  normal  octane,  X^  =  normal 
eicosane,  q°  =  q°'  =  Aq  =  Aq '  =  0,  b°  =  .04,  b°  =  0, 
b°  =  .96,  b°  =  0,  Abx  =  0,  Ab2  =  .96,  Ab3  =  0,  ^  =  .04. 
According  to  Ref.  24,’  the  vapor  pressures  at  40°C  are 
given  by  p°  =  217549.9  mm  Hg,  p2  =  279.4409  mm  Hg,  p°  = 
31.10092  mm  Hg,  p^  =  .0000257  mm  Hg.  Equation  (33)  shows 
that  kx  =  1812.916,  k2  =  2.32867,  k3  =  .259174,  k4  =  2.1«10_ 
Note  that  k^  is  so  large  and  k^  is  so  small  that  methane 
and  normal  eicosane  virtually  play  the  role  of  impermeable 
species . 

Theorem  1  shows  that  S(t)  has  a  unique  solution  for  all 
t  £  0.  Figure  2  shows  the  computed  values  of  x  and  x' 

Ref.  24,  p.  336. 


2.5 
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u  5.0  10.0  15.0 

t 

Fig. 2 — Computed  sizes  of  the  liquid  and  gas  phases  formed  when  0.04  moles  of  methane, 
0.96  moles  of  normal  octane,  0.96  t  moles  of  normal  hexane  and  0. 04 1  moles 
of  normal  eicosane  are  equilibrated  at  40°C  and  120  mm  Hg 
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(the  size  of  the  liquid  and  gas  phases  in  total  number  of 
moles)  plotted  against  t.  Note  that  the  complex  behavior 
of  x  shown  in  Fig.  2  does  not  represent  an  artifact  such 
as  computational  round-off  error.  Intuitively,  one  would 
not  anticipate  this  complicated  double  reversal  in  the 
size  of  the  liquid  phase  as  components  are  added  to  such 
a  simple  system,  involving  no  intraphase  reactions  and 
obeying  the  most  ideal  of  all  possible  laws  for  a  two- 
phase  chemical  system. 

Example  4 

Reference  21  treated  a  primitive  system  with  two 
permeable  species  and  with  an  impermeable  species  repre¬ 
senting  protein,  as  a  tutorial  model  of  human  renal  func¬ 
tion.  In  the  treatment,  Phase  I  represented  plasma  and 
Phase  II  represented  urine.  Far  from  being  simple,  inert, 
and  symmetric,  the  "membrane”  dividing  the  two  phases 
represented  the  aggregated  whole  of  the  complex  secretion 
and  reabsorbtion  processes  involved  in  human  urine  pro¬ 
duction  and  control. 


■nPK-4 
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VII.  CLASSIFICATION  OF  VARIABLE  SYSTEMS 


This  section  discusses  the  behavior  of  the  variable 
system  S(t)  for  t  ^  0  in  terms  of  the  four  possible  types 
of  behavior  for  a  primitive  system  set  forth  in  Theorem  1 
of  Sec.  III. 

Consider  a  variable  primitive  chemical  equilibrium 
system  as  defined  in  Sec.  VI  and  let 

n 

r°  V  ,  °  ,  °  ,  °, 

b  =  lbi  +  q  +q 

i=l 

and 

n 

Ab  =  ^  AIk  +  Aq  +  Aq '  . 

i=l 


In  addition,  we  introduce  the  following  six  quantities: 


°  *  <  n  h° 


to  *  / 


Ah. 


iV‘ 


to  -  Ah 

0 


i  f  q  •  0 


)  TT 

^  i-i  1 


l  f  q  +  Af|  d 

t  f  q°  =  Aq  “  0 

f  Ad  * 

i  I  Ad  -  * 


f  " 

o.  J 

o  -  /  n 


)  k.h, 

/  \  l 

i  *1 


to' 


)  k.A»). 
L  i  i 

i-1 


if  q  '  •  0 

if  q°'  -  0 

if  q"'  +  Aq'  0 

I  [  q"'  =  Aq  '  -II 

i  f  Ad  ■  “ 


Ad'  -  All 

0  i  f  Ad '  = 
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Note:  r  will  not  be  used  if  ct°  =  ®  or  if  Ac  =  Ab;  similarly, 

t'  will  not  be  used  if  a° '  =  ®  or  if  An'  =  Ab . 

The  a’s  are  closely  related  to  Conditions  I  and  II  of 
Sec.  III.  Thus,  Condition  I  for  S(0)  is  equivalent  to  the 
statement  a°  >  b°.  The  r's  represent  certain  critical  values 
of  the  parameter  t. 

Table  1  describes  the  13  possible  relations  between  the 
quantities  a°,  b° ,  cr0*  ,  Ac,  Ab,  and  A o'.  We  assert  that 
every  variable  primitive  chemical  equilibrium  system  satis¬ 
fies  one  and  only  one  of  these  sets  of  four  inequalities. 

Table  1  also  shows  the  types  of  behavior  of  S(t),  as 
t  varies  in  each  of  the  13  cases.  The  phrase  "I  only"  means 
that  S(t)  ’vs  a  unique  solution  and  Phase  II  is  vacuous; 

"II  only"  means  that  S(t)  has  a  unique  solution  and  Phase 
I  is  vacuous.  "Both"  means  that  S(t)  has  a  unique  solu¬ 
tion  and  neither  phase  vanishes.  "Ambig"  (ambiguous)  means 

•k 

that  S(t)  does  not  have  a  unique  solution. 

Finally,  Table  1  indicates,  for  each  of  the  13  cases, 
the  intervals  of  t  in  which  S(t)  has  the  various  types  of 
behavior.  (A  square  bracket  means  that  the  endpoint  in 
question  is  to  be  included;  a  curved  parenthesis  means 

See  Sec.  Ill  for  further  properties  for  each  of  the 
four  types  of  behavior. 
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Table  1 

CLASSIFICATION  OF  VARIABLE  PRIMITIVE  SYSTEMS 


Behavior 


Ac’ 

Ah 

Both 

0  ,®) 

Ac' 

■  Ah 

Both 

'0,r) 

11  only 
"  T  ,  ®  ) 

Ac’ 

•  Ah 

Both 

'0,r') 

I  onlv 

T*  ,®) 

Ac' 

■  Ah 

II  on  1  y 

0  ,  <* ) 

A o' 

■  ■  Ah 

II  omy 
"0,t‘ 

Both 

(T,“) 

Ac' 

Ab 

II  onlv 

'0,T' 

Both 
( T  ,  T  1  ) 

1  onlv 
>’  .®) 

Ac ' 

••  Ah 

I  only 
r0 ,  ®) 

Ac' 

>  Ah 

I  only 
"0,T’  ‘ 

Both 

(r\®) 

Ac' 

-  Ab 

I  only 
:  0  ,  t  1  % 

Both 
(T* , T ) 

11  only 
T  ,  ®  ) 

Ac' 

Ab 

Ambig 

'  0 ,  ®) 


Ac' 

■  Ab 

Ambig 

11 

only 

0  • 

(0, 

“) 

Ac' 

c  Ab 

Amb  ig 

I  only 

■0} 

(0, 

*) 

1J 


Act  Ah 


Ac 1  -  Ah 


Ambig  Both 
: 0 }  (0,®) 
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that  it  is  to  be  excluded.)  Thus,  according  to  Table  1, 
if  a°'  £  b°  <  ct°  and  An  <  Ab  <  A or '  ,  then  S(t)  has  a  unique 
solution  for  all  t  ^  0  and  Phase  I  is  vacuous  for  t  £  t' , 
Phase  II  is  vacuous  for  t  ^  t,  and  neither  phase  is  vacuous 
for  t  <  t  <  T '  . 

Many  readers  may  find  it  instructive  to  derive  Table 
1  from  Theorem  1.  We  omit  this  rather  long  derivation. 
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VIII.  MONOTONICITY  THEOREM 


In  this  section  we  show  that  some  of  the  quantities 
associated  with  a  variable  system  vary  monotonically  with 
t.  This  result  is  of  interest  because,  as  was  illustrated 
in  Example  3  of  Sec.  VI,  other  significant  quantities  can 
behave  very  non-monotonically . 

If  S(t)  has  a  unique  solution,  then  let  x(t) ,  x'(t), 
x^(t),  x^(t),  etc.,  denote  the  values  of  x,  x' ,  x^,  x^ , 
etc.,  for  the  unique  solution  of  S(t). 

THEOREM  2 

If  S(0)  has  a  unique  solution,  then  S(t)  has  a 

x*  (t) 

unique  solution  for  all  t  ^  0,  and  the  ratio  v  1 

x(t) 

(admitting  ®  as  a  possible  value)  is  monotonic 

in  t. 

Proof  of  Theorem  2: 

That  S(t)  has  a  unique  solution  for  all  t  ^  0  follows 
from  Table  1.  It  also  follows  from  Table  1  that  the  set 

T  of  t  >  0  for  which  both  x(t)  >  0  and  x* (t)  >  0  is  an 

interval  of  the  form  (0,®),  (0,t),  (0,t'),  (t,°°),  (t,t'), 
(t',cd),  (t’,t),  or  (0,®).  Furthermore,  it  follows  from 
Table  1  and  from  the  fact  [111  that  x(t)  and  x' (t)  are 
continuous  functions  of  t  that  it  is  sufficient  to  prove 

that  x-  is  monotonic  in  the  interval  T. 

x(t) 
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v  *  (  t  ^ 

But  in  T,  — is  finite- valued , 
x(t) 


and  continuous. 


Thus,  if  this  ratio  were  not  monotonic  in  t,  there  would 
have  to  be  a  real  number  r  and  t^,  t^  (  T  with  t^  <  t^ 
such  that 


x'(t1)  x 1 ( 1 2 ) 

- -  _  _ -  =  r  . 

x(t1)  x ( 1 2 ) 


Moreover,  we  can  choose  t^ 
not  hold  for  any  t  in  (t^, 


x»  (t) 

and  1 9  such  that  — ^  =  r  does 

x(t) 

ty) .  Let  t  =  %(t^  +  t2)  and  let 


Xf  -  h  +  xi^t2^ '  f°r  i=1»  •  •  •  >n 


Xi  =  ^  V xi(ti>  +  xi(t2>)  for  i=l , .  .  .  ,n 


q  =  h  (q(ti)  +  q(t9) 


q'  =  h{  q' (tj^)  +  q ' (t9) 


We  will  show  that  x^,...,x  ,q,  x| ,q  ’  is  an  equi 
librium  composition  for  S(t). 

To  do  this,  observe  first  that 


\ 
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x .  +  x ! 
i  i 


xi(t1)  +  xi(t2), 


x'i(t1)  +  x^(t2) 


=  12  ^(tj)  +  X[(t1^  +  h  t ( 1 2 >  +  x[(t2^) 

=  %  (VV  +  bi(t2j) 


=  b.(t) 


and  that  q  =  q(t)  and  q'  =  q'(t).  Thus,  the  laws  of  mass 
balance  for  S(t)  are  satisfied. 

We  also  have 


x '  = 


x'  (t^  +  x'  (t2) 


=  -  (^(t^  +  x(t2) 


=  rx  . 


It  is  also  clear  that  x  >  0  and  x'  >  0.  Hence  we  may 
write: 
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x!  x  x!  x  x!  (t, )  +  x! (t0) 

1 _ 1  LV  1  1  2 

xi  x'  X.  X*  xi(t1)  +  x^ ( 1 2 ) 

x  x,(t1)x^(t1)  +  x'  (t2)xj^(t2) 
x'  x(t1)xi(t1)  +  x(t2)xi(t2) 

x  rx(t1)kixi(t1)  -I-  rx(t2)kixi(t2) 
x'  x(t1)x1(t1)  +  x(t2)xi(t2) 

=  r  k  =  k 

- ,  i  l 


Thus,  the  laws  of  mass  action  are  satisfied. 

Therefore,  x, ,...,x  ,q,  x' , .  .  . ,x' ,q '  is  an  equilibrium 
1  n  1  n  n 

composition  and,  hence,  the  unique  solution  to  S(t).  Thera- 

.  .  .  .  ^  ^  (t)  x  * 

fore,  x(t)  =  x  and  x'(t)  =  x  so  that  — =  —  =  r,  a  con- 

x(t)  x 


tradiction.  QED. 

COROLLARY  1 

x[(t)  xi(t)  x(t) 

The  functions  - ,  - ,  3 -  and  their  re 

xL(t)  ( t)  b (t) 

ciprocals  are  also  monotonic  in  t. 

Proof  of  Corollary  1: 

From  Eq.  (9)  of  Sec.  I, 


x!(t) 

x^t) 


x'(t) 

=  k.  - 

L  x(t) 
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From  Eq .  (10)  of  Sec.  II, 

xi(t)  x(t)  1 

bt(t)  x(t)  +  k^x' (t)  l  +  k>  x1 (t) 

1  x(t) 

From  Eq.  (6)  of  Sec.  I, 

iM „  i 

b(t)  1  + 

x(t) 


The  reciprocal  of  a  monotonic  function  is  monotonic. 
COROLLARY  2 

If  for  some  i,  =  0,  then  for  that  i,  x 

xl^(t)  are  monotonic  functions  of  t. 

Proof  of  Corollary  2: 

Because 


xt(t)  = 


bT  +  Ab.t 

l  i 


o 


1  +  k 


*'(o 


.x(t) 


1  +  k. 


*'(t) 

Lx(t)  J 


(see  Eq.  (10)  in  Sec.  II),  and  xl^(t)  =  b°  -  x^(t) 


£ED. 


(t)  and 
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COROLLARY  3 

Either  each  of  the  functions  x^(t),  i=l,...,n, 
and  x(t)  is  monotonically  increasing,  _or  each  of  the 
functions  xl^(t)  ,  i=l,...,n,  and  x' (t)  is  monotonically 
increas ing . 

Proof  of  Corollary  3: 

x  *  ( t ) 

If  — is  monotonically  decreasing,  then  we  observe 
x(t) 

that 


\(t)  = 


b.(t) 


1  +  k. 


iKtl 

x(t) 


The  ratio  of  a  monotonically  increasing,  non-negative  func¬ 


tion  to  a  monotonically  decreasing  positive  function  is 


n 


mono 


tonically  increasing.  Since  x(t)  =  q(t)  +  F  x.(t), 

i=l  1 

x(t)  is  also  monotonically  increasing. 

x'(t)  x(t) 

If  — 1 L  is  monotonically  increasing,  then  ■■ v  is 

x(t)  x'(t) 

monotonically  decreasing  and  the  result  follows  by 


reversing  the  roles  of  Phase  I  and  Phase  II. 

Note  the  illustration  of  this  Corollary  in  Fig.  2. 
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IX.  LINEAR  RESPONSES 

The  examples  of  Sec.  VI  indicate  that  the  are 
rarely  linear  functions  of  t.  Under  what  circumstances 
will  linearity  hold? 

THEOREM  3 

Assume  that  S(0)  has  a  unique  solution  with 
neither  phase  vacuous;  then  x^(t),  i=l,...,n,  will 
all  be  linear  functions  of  t,  0  ^  t  <  ®,  precisely 
if  either: 

Case  A:  There  are  constants  u,  u'  such  that 

=  ux^(O)  +  u’x|(0)  for  1=1,..., n  ,  (40) 

=  uq°  ,  (41) 

V  =  u' q° 1  ,  (42) 

in  which  case: 


xi(t)  -  (1  +  tu)x1(0)  ,  (43) 

x!(t)  =  (1  +  tu')x!(0)  ,  (44) 
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x(t)  =  (1  +  tu)x(0)  , 

(45) 

x'  (t)  =  (1  +  tu)x' (0)  , 

(46) 

x.(t)  =  x.(0)  , 

(47) 

x!(t)  =  x!(0)  ; 

(48) 

or  , 


Case  B:  The  constant  r 


*'(0) 

x(0) 


is  such  that 


l 


i=l 


Ab.  (1-k.) 
i  1 


1  +  rk. 

l 


Aq 


5 


(49) 


in  which  case 


b.(t) 

xi(t)  =  1  +  rk.  ’ 


rk.b. (t) 
L 


sim  =  r 

i(t) 


(50) 


(51) 


(52) 


-  —  v*  »  ■r» 


m. 
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Proof  of  Theorem  3: 

Since  S(0)  has  a  unique  solution,  so  does  S(t)  for 
t  ^  0.  Assume  first  that  each  x^(t)  is  a  linear  function 
of  t.  Since  b^t)  ,  q(t),  q'(t),  b(t)  are  all  by  defini¬ 
tion  linear  functions  of  t,  then 

x[(t)  =  bt(t)  -  xt(t)  , 
n 

x(t)  =  £  xt(t)  +  q(t)  , 

i=l 

x' (t)  =  b(t)  -  x(t) 

are  also  linear.  By  hypothesis,  x(0)  >  0;  hence,  for  any 
t  *  0, 


(0)  +  x(2tjj  *  *121  >  o  . 

Similarly,  x* (t)  >  0  for  t  s  0.  We  can  then  assert  the 
law  of  mass  action: 


x[(t)  kixi(t) 

-  - - for  n=l , .  .  .  ,n  and  t  ^  0  . 

x’(t)  x(t) 


(53) 
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Hence , 


x^(t)x(t)  =  k^x^(t)x'(t)  for  n-l,...,n  and  t  2  0 

(54) 

A  linear  function  defined  for  t  ^  0  can  be  extended 
linearly  to  the  whole  real  line.  If  this  is  done  for 
x^(t) ,  x|(t),  x(t) ,  x' (t)  ,  then  Eq.  (54)  remains  valid  for 
all  real  t,  and  represents  two  ways  of  factoring  the  same 
quadratic  function  of  t  into  linear  factors.  But  since 
the  linear  factors  correspond  to  the  roots  of  the  quadratic, 
this  factorization  is  unique  except  for  multiplicative  con¬ 
stants  and  the  order  of  the  factors.  Thus  we  have  shown 
that  either  x(t)  is  proportional  to  x^(t)  or  x(t)  is  pro¬ 
portional  to  x'(t).  Since  x(t)  >  0  for  t  >  0,  we  have 

shown  that:  either  for  every  i 

xi(t)  xi(0) 

■  —  is  a  constant,  viz.,  -  ;  (55) 

x(t)  x(0) 

or 


him 

x(t) 


is  a  constant,  viz., 


him  =  r 

x(0) 


(56) 


In  the  first  case,  we  have  from  Eq.  (53)  that  for 


every  i, 


x[(t)  x’(0) 

-  is  a  constant,  - 

x'(t)  x’ (0) 

Define 


d_  x(t) 
dt  x(0) 


and 


.  d  x'(t) 
a'  - - ; 

dt  x'(0) 

since  x  and  x'  are  linear,  u  and  u*  are  constants. 


£b 


i  "It  bi(t)  =  +  xi(b>) 


2^-  +x!(0)  ^Vux.(O) 


m  f-  x.(0) 

d  ^  1  x(0)  1  x'(0) 


Then 


+  u'x 


which  is  Eq.  (40) 


Also , 
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Aq  -  q(t) 


u 


> 


which  is  Eq.  (41).  The  proof  of  Eq.  (42)  is  similar. 

Then  Eqs.  (43)  through  (48)  follow  by  simple  manipulations. 
In  the  second  case,  Eq.  (56),  we  have 


x(t) 


r+1 

r+1 


x(t) 


x(t)  +  x1  (t) 
r+1 


_  EsI 

r+1 


n 

q (t)  +  q' (t)  +  I  b. (t) 

_ i=l 

r+1 


(58) 


Now  x(t)  must  be  the  root  of  the  characteristic  function 
of  S(t) ;  that  is , 
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q(0  " 

I -  +  )  7 - -  -  1  . 

x(t)  1=1  x(t)  +  k^x' (t) 

Multiplying  both  sides  of  this  by  x(t)  gives 
"  bi(t)  . 

<i(t)  +  1  ritir  =  x(c)  • 

.  _  i 

i=l 

Replacing  x(t)  by  its  value  from  Eq.  (58),  differ¬ 
entiating  with  respect  to  t,  and  simplifying,  gives  Eq. 
(49);  then  Eqs.  (50),  (51),  and  (52)  follow  easily. 

Finally,  it  is  an  algebraic  exercise  to  show  that,  if  the 
condition  for  either  case  A  or  case  B  holds,  then  the  x^(t) 
have  the  claimed  behavior  and  hence  are  linear.  QED. 

Note  that  Case  A  of  the  theorem  corresponds  simply  to 
expanding  the  two  phases  separately,  without  changing  the 
mole  fractions.  One  type  of  variable  system  which  satis¬ 
fies  Case  B  is  given  by  Aq  =  Aq'  =0  and  =  0  for  all 
i  for  which  k^  t  1.  That  is,  if  a  system  is  varied  by 
changing  only  the  inputs  of  "undriven"  substances  (i.e., 
substances  for  which  =  1) ,  then  the  response  is  linear. 
Here  the  x^(t)  and  x^(t)  are  constant  for  the  substances 
whose  inputs  are  not  altered. 
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For  systems  in  which  Aq  =  Aq*  =0  and  for  which 
b°  >  0,  for  all  i,  we  may  simplify  Eq.  (49).  Observe 
from  Eq.  (10)  that 


1  +  rkt 


xt(0) 


J 


hence , 


1  -  k. 

_ i_ 

1  +  rk^ 


x1(0) 

~b?“ 

1 


(1  -  kt) 


xt(0)  -  x[(0) 


so  Eq.  (49)  becomes 


n 


l 


=  0  . 


(59) 


W1 
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